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ABSTRACT 



A method for determining a suggested value for a cancel- 
lation parameter for a dual energy decomposition includes 
obtaining a first energy level image of internal structure, 
obtaining a second, lower, energy level image of the internal 
structure, and iteratively processing the images to determine 
a provisional value for the cancellation parameter. The 
iteration includes varying a cancellation parameter in a 
cancellation equation, obtaining a structure cancelled image 
from the first and second energy level images according to 
the cancellation equation, and evaluating a cancellation 
metric from the structure cancelled image. The provisional 
cancellation parameter may then be chosen (e.g., as the 
value that approximately minimizes a variance cancellation 
metric). Further iterations may be performed around the 
provisional cancellation parameter to refine the provisional 
cancellation parameter into a final cancellation parameter. 

36 Claims, 2 Drawing Sheets 
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DUAL ENERGY DECOMPOSITION USING BRIEF SUMMARY OF THE INVENTION 

AUTOMATICALLY DETERMINED . . . . . . , ... , . - A 

CANCELLATION PARAMETERS Apreferred embodiment of the present invention provides 

a method tor determining a suggested value tor a cancella- 

BACKGROUND OF THE INVENTION 5 lion parameter for a dual energy decomposition. The method 

Tne present invention relates to medical diagnostic X-ray *f lu <| cs °"f™?B a ™W level image of internal 

imaging. In particular, the present invention relates to dual J ™ re ' ° bt f m,ng * ' enef8y kve ' lma 8 e of 

energy decomposition for tissue specific imaging using a internal structure, and iterativcly processing the unages 

* - \ j j * *• * u - * u* • ii to determine a provisional value for the cancellation param- 

computer assisted detection technique to obtain a caacella- T ^- i , . 

lion parameter 10 e particular, the iteration includes varying a cancella- 

m , . " , . . , , uon parameter in a cancellation equation, obtaining a struc- 

Today, doctors and technicians commonly have access to , ure cmcdkd j &om , he ^ and ene ^ 

very sophisticated medical diagnostic X-ray imaging i ma ges according to the canceUauon equation, and evaluat- 

devices. Typically dating the operation of an X-ray imaging mg a metric from me stiucture ^ ; 

device an X-ray source emits X-ray photons under very ^ provisional cance llation parameter may then be chosen 

controlled circumstances. The X-ray photons travel through « ( as ^ value tha , approximatel y mininliz e S a variance 

a region of interest (ROI) of a patient under examination and CMlce iiation metric). Further iterations may be performed 

impinge upon a detector. In the past X-ray imaging devices around ^ provisional canceU ation parameter to refine the 

employed rudimentary film based detectors. However, provisional cancc ii ati on parameter into a final cancellation 

recent developments have led to solid state detectors com- parameter 

prised of a grid of discrete detector elements that individu- 20 „ * it . , 

ally respond to exposure by X-ray photons. Regardless of the P*** invention may be embodied in a 

the detector used, however, the goal remain s the same, me f'= al ^ osi ^ "nagmg processing system. The system 

namely, to produce a clear resultant image of preselected mcludes a P rocessi ?g circ ^ • memory coupled to the 

structures of interest (e.g., specific types of tissues) within pro«^g circuit for storing a first energy level image of 

the ROI M internal structure and a second, lower, energy level image of 

_ ' . . . . . ■ the internal structure. The memory stores instructions for 

Tnere is an inherent difficulty associated with producing execution b ^ ^ tQ a lish ^ st no , ed 

a clear resultant image ftowever. In particular, because the above _ fa articul me 

instructions iterative ly vary a can- 

X-ray photons travel through the entire patient^ the image ceUation eter m , cancellation equation> obtain a 

formed on the detector is asuperposition of all the anatomic canceJled ^ from the flrst and secBaA 

structures through which X-ray photons pass, including the ^ j iccordin to the canc6U ation equation, and 

preselected structures of interest. The superposition of ana- & cancellation metric from lhe structure canceUed 

tomic structures is sometimes referred to as anatomic ■ • * ** *u i * • • 1 « 

r . • . • rmage. The instructions then select a provisional cancella- 

noise . The effect of anatomic noise on the resultant image , ■ * nn „„ t „ u nr . a A -n«*-« ~ ♦ ■ 

, , . , , t , _ & Uon parameter based on the cancellation metric, 
is to produce clutter, shadowing, and other obscuring effects 

that render the resultant image much less intelligible than the 35 BRIEF DESCRIPTION OF THE DRAWINGS 
ideal clear resultant image. 

Past attempts to reduce the effects of anatomic noise . FI 9* 1 fllustrates a hi g h level diagram of an X-ray 

included, for example, "dual-energy" imaging. When imaging system. 

employing dual-energy imaging, a doctor or technician 4Q FIG. 2 illustrates a flow diagram for the operation of the 

acquired an image at high average X-ray photon energy, and X-ray imaging system, 
an image at low average X-ray photon energy. Because 

different internal structures absorb different X-ray photon DETAILED DESCRIPTION OF THE 

energies to different extents, it was possible to combine the INVENTION 

two resultant images to suppress anatomic noise, according <5 Xumiag now to FIG x that figure ilhlstrates arj x . ray 

t0: imaging system 100. The imaging system 100 includes an 

X-ray source 102 and a collimator 104, which subject 

SB(x,y)-ex P [log (H(x,y))-w logfL^y))],^ <i), structure under examination 106 to X-ray photons. As 

examples, the X-ray source 102 may be an X-ray tube, and 
where SB is the decomposed image achieved through the so the structure under examination 106 may be a human 
log subtraction at a specific cancellation parameter w, patient, test phantom or other inanimate object under test. 
H(x,y) is an image obtained at high energy, and L(x,y) The X-ray imaging system 100 also includes an image 
is an image obtained at low energy. By varying w, SB sensor 108 coupled to a processing circuit 110. The pro- 
becomes a decomposed image of either soft tissue or of cessing circuit 110 (e.g., a microcontroller, microprocessor, 
bone. 55 custom ASIC, or the like) couples to a memory 112 and a 
However, in the past, users of the previously mentioned display 114. The memory 112 (e.g., including one or more 
decomposition technique had to vary the cancellation of a hard disk, floppy disk, CDROM, EPROM, and the like) 
parameter, w, manually through trial and error. The resulting stores a high energy level image 116 (e.g., an image read out 
manual variation of the cancellation parameter was time from the image sensor 108 after 110-140 kVp 5 mAs 
consuming and hindered the workflow in the clinical envi- 60 exposure) and a low energy level image 118 (e.g., an image 
ronment. Furthermore, the final value chosen for the can- read out after 70 kVp 25 mAs exposure). The memory 112 
cellation parameter was not always the one that provided the also stores instructions for execution by the processing 
best cancellation of bone or soft tissue given the high and circuit 110, as explained below, to cancel certain types of 
low energy images. structure in the images 116-118 (e.g., bone or tissue 
A need has long existed in the industry for a method and 65 structure). A structure cancelled image 120 is thereby pro- 
apparatus for dual energy decomposition that addresses the duced for display. The image sensor 108 may be a flat panel 
problems noted above and previously experienced. solid state image sensor, for example, although conventional 
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film images stored in digital form in the memory 112 may 
also be processed as disclosed below as well. 

The processing circuit 110 uses a log subtraction tech- 
nique to decompose spatially registered images from the 
structure under examination 106 into constituent materials 
(e.g., bone and soft tissue in chest X-ray images). To this 
end, a high energy level image and a Low energy level image 
are obtained and processed according to: 

SB(x,y)- & xp[log(H(x,y))-w log(L(X x,y))],(0<w<l) (1) 

where SB is the decomposed image achieved through the 
log subtraction at a specific cancellation parameter w, 
H(x,y) is the high energy level image, and L(x,y) is the 
low energy level image. By varying w, SB becomes a 
decomposed image of either soft tissue or of bone. 
The processing circuit 110 iteratively varies the cancel- 
lation parameter w in Equation (1), obtains a structure 
cancelled image (i.e., SB(x,y)) according to Equation (1), 
and evaluates a cancellation metric from the structure can- 
celled image. The processing circuit then selects a provi- 
sional cancellation parameter based on the cancellation 
metric. 

More specifically, the processing circuit 110 minifies 
H(x,y) and L(x,y) (i.e., reduces their size through neighbor- 
hood averaging, for example), and uses Equation (1) at 
varying w to create multiple images of SB(x,y). The param- 
eters for Equation (1) may be selected based on empirical 
knowledge of image acquisition parameters. For example, if 
H(x,y) and L(x,y) are obtained at 140 kVp, 5 mAs and 70 
kVp, 25 mAs, a starting cancellation parameter, w, may be 
0.2 and the step size may be 0.05. On the other hand, if 
H(x,y) and L(x,y) are obtained at 110 kVp, 5 mAs, and 70 
kVp, 25 mAs, a starting value of the cancellation parameter 
may be 0.4 and the step size may be 0.04, The cancellation 
parameter may then be varied between a range of 0.3 and 
0.95, or other ranges empirically determined. 

Note that no particular magnitude of difference in the first 
energy level and the second energy level is required. The 
iterative determination of the cancellation parameter pro- 
ceeds independently of the energy levels used. 

The images H(x,y) and L(x,y) may vary depending on the 
type of structure sought to be cancelled. For example, if 
bone is to be cancelled, then H(x,y) and L(x,y) may be 
64x64 pixel images of a portion of the spine (or other 
predominately bone structure). If tissue is to be cancelled, 
then H(x,y) and L(x,y) may be 256x256 images of a lung 
region between ribs (or other predominately tissue area). 

At each step, a cancellation metric (e.g., variance iD 
SB(x,y)) is determined. The variance will be lowest at the 
value of w achieves the most cancellation of bone (for 
H(x,y) and L(x,y) images of bone structure) or of tissue (for 
H(x,y) and L( x >y) images of tissue structure). The lowest 
variance is selected as a provisional cancellation parameter 
and may be used to cancel similar structure in other images. 
Preferably, however, the processing circuit 110 performs one 
or more additional iterations at smaller step sizes around the 
provisional cancellation parameter. 

For example, the new step size may be 0,01, and new 
SB(x,y) images are then computed according to Equation (1) 
on either side of the provisional cancellation parameter. The 
processing circuit 110 determines the variance in each 
SB(x,y) to select a final cancellation parameter that best 
cancels bone or tissue. It is noted that higher resolution 
H(x,y) and L(x,y) images may be used during the additional 
iterations. 

The process identified above may further automated 
through the use of a histogram technique. In particular, the 
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processing circuit 110 may employ a pattern recognition 
algorithm to identify the chest cage, spine, and ribs. In the 
central region of a standard posterior-anterior image, a larger 
X-ray absorption region in the middle of horizontal rib lines 

5 (or the average of horizontal rib lines) represents the spine. 
Away from the spine, periodic high transmission intensities 
in vertical lines represent unobstructed lung, while a lower 
X-ray transmission represents ribs. 

Thus, an intensity histogram may be constructed to iden- 

10 tify two lungs based on pixel intensity. A profile drawn 
across the lungs will include two peaks (one for each lung), 
and a valley (corresponding to the spine). The peaks and 
valleys then identify structure from which H(x,y) and L(x,y) 
may be taken as predominately bone (e.g., the spine) or 

35 tissue (e.g., the lung). 

It is also noted that the processing circuit 110 may first 
remove background radiation from the H(x,y) and L(x,y) 
images before selecting a cancellation parameter. To this 
end, the processing circuit 110 may use regions of air (e.g., 

20 regions between the chest and arm, or above the shoulders) 
identified, for example, in an image or in the intensity 
histogram to determine a background radiation level. The 
background radiation level may then be subtracted from 
H(x,y) and L(x,y). 

25 Turning next to FIG. 2, that figure present a flow diagram 
200 that shows the processing steps that occur at the 
processing circuit 110. First, the processing circuit 110 
acquires an image (202) an prepares an intensity histogram 
(204). Bone and tissue structure, as well as background 

30 radiation are identified (steps 206-208). Subsequently, the 
processing circuit 110 obtains H(x,y) and L(x,y) through a 
high and low energy imaging selection (step 210) and 
removes background radiation, if desired, from H(x,y) and 
L(x,y) (step 212). 

35 Continuing at step 214, the processing circuit iteratively 
varies the cancellation parameter, w, according to a first step 
size and evaluates Equation (1) to obtain a structure can- 
celled image SB(x,y). The variance in SB(x,y) is determined 
(step 216) and a provisional value for the cancellation 

40 parameter is selected (step 218). Subsequently, at step 220, 
the processing circuit 110 reduces the step size, and itera- 
tively evaluates Equation (1) by varying the cancellation 
parameter around the provisional value (step 222) by the 
reduced step size. 

45 Similarly, at step 224, the processing circuit determines 
the variance in SB(x,y). At the end of the iterations, the 
provisional cancellation parameter may be updated (step 
226) and additional iterations performed with smaller step 
sizes. Alternatively, the processing circuit may instead select 

50 a final value for the cancellation parameter that yields the 
minimum variance in SB(x,y) (step 228). 

Thus, the present invention provides an automated tech- 
nique for selecting a value for a cancellation parameter that 
best eliminates structure (e.g., bone or tissue) from an 

55 image. The final value of the cancellation parameter may 
thus be obtained without significant trial and error delays, or 
subject to possible operator biases. The cancellation param- 
eter may then be used to cancel the same structure in 
subsequently obtained images. The software underlying the 

60 processing steps shown in FIG. 2 for the X-ray imaging 
system 100 may be stored on a computer readable medium 
(such as the memory types identified above) for convenient 
distribution. 

While the invention has been described with reference to 
65 a preferred embodiment, those skilled in the art will under- 
stand that various changes may be made and equivalents 
may be substituted without departing from the scope of the 
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invention. In addition, many modifications may be made to 14. An image processing system for medical diagnostic 

adapt a particular step, structure, or material to the teachings imaging, the image processing system comprising: 

of the invention without departing from its scope. Therefore, a processing circuit; 

it is intended that the invention not be limited to the a memory coupled to the processing circuit, the memory 

particular embodiment disclosed, but that the invention will S comprising storage for a first energy level image of 

include all embodiments falling within the scope of the internal structure and a second energy level image of 

appended claims. the internal structure at an energy level lower than the 

What is claimed is: first energy level image, 

1. A method for determining a cancellation parameter for the memory further comprising instructions for execution 
a dual energy decomposition, the method comprising: 10 by t& e processor for: 

obtaining a first energy level image of internal structure; iteratively varying a cancellation parameter in a can- 

4 . . , , . j. . . , cellation equation, obtaining a structure cancelled 

obtaining a second energy level image of the mternal image from the first ^ energy [cvcl images 

structure at an energy level lower than the first energy according to the cancellation equation, and evaluat- 

level image; ^ j ng a cancellation metric from the structure cancelled 

iteratively varying a cancellation parameter in a cancel- image; and 

lation equation, obtaining a structure cancelled image selecting a provisional cancellation parameter based on 

from the first and second energy level images according the cancellation metric. 

to the cancellation equation, and evaluating a cancel- 15. The image processing system of claim 14, wherein the 

lation metric from the structure cancelled image; and 2Q memory further comprises instructions for obtaining the first 

selecting a provisional cancellation parameter based on and second energy level images from an X-ray detector, 

the cancellation metric. 16. The image processing system of claim 14, wherein the 

2. The method of claim 1, wherein the internal structure internal structure is predominately bone structure, and the 
is predominately bone structure, and the cancellation metric cancellation metric is a bone cancellation metric. 

is a bone cancellation metric. 25 17 - ^ ^ag 6 processing system of claim 16, wherein the 

3. The method of claim 2, wherein the internal structure internal structure comprises a portion of a spinal column, 
comprises a portion of a spinal column. 18 - The ™age processing system of claim 14, wherein the 

4. The method of claim 1, wherein the internal structure internal structure is predominately tissue structure, and the 
is predominately soft tissue structure, and the cancellation cancellation metric is a tissue cancellation metric, 
metric is a soft tissue cancellation metric. 30 19 - T^ method of claim 18, wherein the internal structure 

5. The method of claim 4, wherein the internal structure comprises a portion of a lung. 

comprises a portion of a lung. 20 - ^ ima e e processing system of claim 14, wherein the 

6. The method of claim 1, wherein the cancellation cancellation equation is: 
equation is: 

SB(x,y)=exp[logO(H(x,y))-w log(L(x,y))], 

35 

SB(x,y)=exp[log(H(x,y))-w logfLfoy))!, and w is the cancellation parameter. 

21. The image processing system of claim 14, wherein the 

and w is the cancellation parameter. cancellation metric is variance in the structure cancelled 

7. The method of claim 1, wherein the first and second image. 

energy level differ by at least 40 kVp. 22 ^ image processing syslem of claim 14) wherein the 

8. The method of claim 7, wherein the first and second instructions for iteratively varying comprise instructions for 
energy level differ by at least 70 kVp. iteratively varying the cancellation parameter by a first step 

9. The method of claim 1, wherein the cancellation metric sizc over a predetermined range. 

is variance in the structure cancelled image. 23. The image processing system of claim 22, wherein the 

10. The method of claim 1, wherein iteratively varying memory stores instructions for: 

comprises iteratively varying the cancellation parameter by sekcti a second ^ smaHer than ^ ^ gi 

a first step size over a predetermined range. 4 . , . lU ^ A . iL 

^ ™ r ° .. , iterative lv varvins the cancellation narameter in the can - 

11. ine meinoa or claim iO, runner comprising, suose- 1t * ' r , . "... „ . 

quent to the selecting step, the steps of: cellatl0 ° ^T" "Tl . pr0V ^° nal canc ( ella lon 

^ ° r r parameter by the second step size, obtaining a structure 

selecting a second step size smaller than the first step size; 5Q cancelled image from the first and second energy level 

iteratively varying the cancellation parameter in the can- images according to the cancellation equation, and 

cellation equation around the provisional cancellation evaluating the cancellation metric from the structure 

parameter by the second step size, obtaining a structure cancelled image; and 

cancelled image from the first and second energy level selecting a final cancellation parameter based on the 

images according to the cancellation equation, and 55 cancellation metric. 

evaluating the cancellation metric from the structure 24. The image processing system of claim 23, wherein the 
cancelled image; and memory further stores instructions for subsequently cancel- 
selecting a final cancellation parameter based on the ing structure in additional images using the cancellation 
cancellation metric. equation and the final cancellation parameter. 

12. The method of claim 11, further comprising the step 60 25. The image processing system of claim 14, wherein the 
of subsequently canceling structure in additional images memory further stores instructions for subsequently cancel- 
using the cancellation equation and the final cancellation ing structure in additional images using the cancellation 
parameter. equation and the provisional cancellation parameter. 

13. The method of claim 1, further comprising the step of 26. The image processing system of claim 14, wherein the 
subsequently canceling structure in additional images using 65 first and second energy level differ by at least 40 kVp. 
the cancellation equation and the provisional cancellation 27. The image processing system of claim 26, wherein the 
parameter. first and second energy level differ by at least 70 kVp. 
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28. A computer program product for structure cancellation 
in images of internal structure, the product comprising: 

a storage medium readable by a processing circuit and 
storing instructions for execution by the processing 
circuit for: 5 
iteratively varying a cancellation parameter in a can- 
cellation equation, obtaining a structure cancelled 
image according to the cancellation equation from a 
first energy level image obtained at a first energy 
level and a second energy level image obtained at a 10 
second energy level different than the first energy 
level, and evaluating a cancellation metric from the 
structure cancelled image; and 
selecting a provisional cancellation parameter based on 
the cancellation metric. 15 

29. The computer program product of claim 28, wherein 
the internal structure is predominately one of tissue and bone 
structure, and the cancellation metric is variance in the 
structure cancelled image. 

30. The computer program product of claim 28, wherein 20 
the cancellation equation is: 

SB(x,y)=e X p[log(H(x,y))-w log(L(x,y))], 

and w is the cancellation parameter. 2 5 

31. The computer program product of claim 28, wherein 
the instructions for iteratively varying comprise instructions 
for iteratively varying the cancellation parameter by a first 
step size over a predetermined range. 
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32. The computer program product of claim 31, wherein 
the storage medium further stores instructions for: 

selecting a second step size smaller than the first step size; 

iteratively varying the cancellation parameter in the can- 
cellation equation around the provisional cancellation 
parameter by the second step size, obtaining a structure 
cancelled image from the first and second energy level 
images according to the cancellation equation, and 
evaluating the cancellation metric from the structure 
cancelled image; and 

selecting a final cancellation parameter based on the 
cancellation metric. 

33. The computer program product of claim 32, wherein 
the storage medium further stores instructions for subse- 
quently canceling structure in additional images using the 
cancellation equation and the final cancellation parameter. 

34. The computer program product of claim 28, wherein 
the storage medium further stores instructions for subse- 
quently canceling structure in additional images using the 
cancellation equation and the provisional cancellation 
parameter. 

35. The computer program product of claim 28, wherein 
the first and second energy level differ by at least 40 kVp. 

36. The computer program product of claim 35, wherein 
the first and second energy level differ by at least 70 kVp. 

***** 
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